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Using an example of a Raman active medium we describe how a common nonlinear process of 
four-wave mixing can be used to induce strong coupling between the spatial and temporal degrees 
of freedom in optical waves. This coupling produces several unexpected effects. Amongst those are 
cascaded excitation of multiply charged optical vortices, spatial focusing in a nonlinear ly defocusing 
medium and generation of helically shaped spatio-temporal optical solitons. 



The nonlinearly mediated interaction between fre- 
quency harmonics and the associated effect of frequency 
conversion have laid the foundation of nonlinear optics 
[H , and have had an important effect in areas such as fluid 
dynamics and plasma physics [2[, and - most recently - 
physics of ultracold matter Four- wave mixing is the 
most common type of nonlinear wave interaction seen 
in both crystals with inversion symmetry and isotropic 
media. In the context of nonlinear optics the four-wave 
mixing process is typically mediated by Kerr and Raman 
nonlinearities [l|. 

Nonlinearity can be used not only to control the spec- 
tral, and hence the temporal, properties of waves, but 
also to transform the linear [l| and orbital angular mo- 
menta of waves 0, 0, 0, HI- In particular, experiments 
with the second harmonic generation by optical vortices 
0] have demonstrated that two photons with frequency 
uj and an orbital angular momentum (OAM) quantum 
number (or vortex charge) / produce a single photon 
with frequency 2uo and OAM quantum number 21 [5[. 
Analogous OAM conversion rules have been reported for 
soliton-like beams within the second harmonic generation 
model [4|, ML for degenerate four- wave mixing in Kerr-like 
materials |8[ , and for three- wave Raman resonant process 
with higher order Bessel-beams [6[ . Though the interplay 
between spatial and temporal dynamics in nonlinear me- 
dia has recently become an active research topic, see, e.g. 
[sj], it has failed so far to exploit many of the promising 
consequences of the strong spatio-temporal coupling seen 
in frequency conversion experiments using beams with 
non-zero OAM. 

In this work we suggest a method for manipulating 
the frequency and angular harmonics of light opening 
new opportunities for simultaneous spatial and temporal 
shaping of optical waves. We will demonstrate that the 
exploitation of the coupling between spatial and tempo- 
ral degrees of freedom can lead to such exciting effects as 
cascaded vortex generation, strong spatial focusing in- 
duced by self-defocusing nonlinearity and the generation 
of spatio-temporal helical beams in the solitonic and non- 
solitonic regimes. 

One of the most efficient nonlinear optical processes is 
the Raman mediated four- wave mixing in gases, which 
easily produces dozens of spectral lines. In particular, it 



has been demonstrated that the two frequency excitation 
of a Raman transition away from the resonance can gen- 
erate, by means of the cascaded four-wave mixing, broad 
frequency combs and induce nonlinearity related chirp 
[ToL fill . EH- The nonlinear chirp can be compensated 
for by the intrinsic group velocity dispersion (GVD) of 
the material resulting in the generation of trains of ultra- 
short (close to single-cycle) pulses [Io|, EH EI| • Note, that 
the generation of the frequency combs with this tech- 
nique relies not on the Raman gain, but on a four- wave 
mixing process which dominates the wings of the Ra- 
man line and induces strong effective Kerr nonlinearity 
fl3| | . The starting idea for the results presented below 
is that the cascaded off-resonant Raman process, as in 
[l(j, is initiated under conditions where one of the two 
(frequency detuned) driving fields is a singly charged vor- 
tex beam. In this case, the phase dependent nonlinear 
coupling between the Raman side-bands (the same cou- 
pling, as that which drives frequency conversion) causes 
cascaded generation of multiply charged vortices. A di- 
mensionless model describing the evolution of the side- 
bands is 0, El 
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where n = -TV, . . . , 0, . . . , N and A = d^+dy. E n are the 
dimensionless amplitudes of the sidebands, such that the 
total field is given by E tot = £) n E n (x, y, z)e innt ~ iKnZ , 
where Q n (cj + ncj mod ) /(j mod . cj mod = ui - cj is 
the modulation frequency (which is the frequency differ- 
ence between the two driving fields) and n = 0,1. The 
physical frequencies and wavenumbers are represented by 
the lower case letters uo n and fc n , whilst their dimen- 
sionless equivalents - by the upper case: Q n and K n . 
The dimensionless time t is measured in units of l/umodi 
the propagation coordinate z is in units of L, and the 
transverse coordinates (x,y) are in units of ^/Lc/uoq. 
K n = {uoQ J rnuo rno d)L/c is the scaled free space wavenum- 
ber. Here, L = (^^cjoA/"^!) -1 characterizes the coupling 
length over which power is transferred between neighbor- 
ing side-bands in the absence of dispersion, r] ~ 376, J\f 
is the density of molecules and b is a coefficient char- 
acterizing the material dependent coupling between the 
sidebands [To[, which weak frequency dependence is ne- 
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FIG. 1: (Color online) The smaller panels show the transverse 
profiles of the sideband phases (aig(E n (x,y))). The larger 
panel shows the x- distribution of intensities, \E n (x, y — 13) | 2 . 
The propagation distance is z = 0.2. The window size in 
each plot is 25 x 25 dimensionless units (which corresponds to 
c± 600 x 600/ira). The other parameters are: ji = 5, j = — 1, 
P2 = 0.005, A/" = 23. The initial conditions are: Aq — y/2, 
Ai = 1, wo = 2y / 5, = a/T0 (see text for more details). 



glected for simplicity. 

L varies from 1 to a few mm for D2 and H2 gases 
[Tot , so that the one unit of x corresponds to a few tens 
of microns. the Raman coherence, is responsible for 
the coupling between the side-bands. In the adiabatic 
approximation [3, EH 
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where /i = |cj m0( i — oj r |/(|6|7o) is the scaled modulus of 
the detuning of the modulation frequency from the Ra- 
man frequency uo r and j = sign(u mo d — uo r ) is the de- 
tuning sign. \q\ varies from to 1/2 for |5| varying from 
to 00. j = +1 or — 1 correspond, respectively, to the 
self-focusing or self-defocusing nonlinearity induced by 
the Raman coherence 13, Il4l [TBI, [lij . E n \fT$ are the 



dimensional amplitudes of the harmonics. For D2 and 
H2 gases /jl = 1 corresponds to Iq ~ 0.1GW/cm 2 , pro- 
vided \uomod — w r \ ~ 1GHz. The frequency dependence of 
the propagation constant j3{uo) can be fitted with a poly- 
nomial, which in the discretized frequency space gives 
/3(co n ) = (3 n = YliLoPi nl - Without any loss of generality 
the po,i coefficients can always be set to zero by phase 
rotations P2 and Vi>2 characterize GVD and higher 
order dispersions respectively. 

In order to derive the rules for the OAM conversion we 
assume that the system is pumped with two neighboring 
sidebands (n = 0,1) having phases 0o and 0i, where 
(j) n = to n t - (K n - (3 n )z + InO, = arg(x + iy) is the 
polar angle and l n is the OAM quantum number of the 
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FIG. 2: (Color online) Intensity distribution of the total field 
\E to t\ 2 after different propagation distances z for t = 0. Az- 
imuthal compression is seen for z = 20 and subsequent recom- 
pression for z = 97. The nonlinearity is defocusing, j = —1. 
wo = 20 \/5, wi = 10 a/T0. The other parameters are the same 
as in Fig.[T] The window size is 250 x 250 dimensionless units. 





FIG. 3: (Color online) Spectral (a) and temporal (b) profiles 
of the total field z = 10 (black lines) and z = 20 (red lines). 
To plot (a,b) we used the cut through the helix taken at the 
point of maximal intensity \Etot\ 2 ■ Other parameters are the 
same as in Fig. [2] (c) Isointensity (x, 2/, £)-plot at 50% of 
maximum for z = 20. 



n's sideband. The conditions for the phase matched ex- 
citation of the nearest Stokes and anti-Stokes sidebands 
are 2</>o — 4>i = 0-i, 20i — <po = 4>2- For an arbitrary 
n these are generalized as 2(j) n = </> n +i + 4> n -i- While 
the ^-dependent part of the above condition can not be 
made zero in dispersive materials (/3 n 7^ 0), the t and 
dependent parts are nulled exactly for uo n = uoq + ncu mo( i 
and l n = Iq + n(li — lo). Therefore, the OAM conversion 
rule is analogous to the one for the frequencies and h—lo 
can be considered as the seeded modulation of OAM. In 
a particular example considered below Iq = and Zi = 1, 
so that l n = n. For sufficiently strong input fields 
Raman four-wave mixing should develop in a cascaded 
manner and excite higher-order Stokes and anti-Stokes 
fields carrying progressively larger OAM quantum num- 
bers. This process is illustrated in Fig. 1, where we show 
the transverse profiles of the phases of the side-bands 

found by numerical integration of Eqs. ([Tj), ([2j) with 47 

2 / 2 

sidebands and the initial conditions = Aoe~ r 
E x = A x re ie - r *K, E n ^ ,i = (r 2 = x 2 + y 2 ). One can 
see that the generated Stokes (n < 0) and anti-Stokes 
(n > 0) components carry progressively increasing vor- 
tex charges, which can be inferred from the number of 
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jumps between the blue (zero phase) and red (phase 2tt) 
colors (see areas around the beam centers in Fig. 1). The 
efficiency of the OAM conversion (in terms of the phase 
matching requirements) is the same as the efficiency of 
the frequency conversion. However, when dealing with 
spatially inhomogeneous beams one should provide suf- 
ficient initial overlap of the two pump fields to build a 
strong coherence q. This is achieved by using a vortex 
free Gaussian beam wider than the one with a vortex 
(w > Wi). 

The above result suggest that the methods of frequency 
domain wave synthesis El, EH have the potential to 
be extended to the angular momentum components, and 
so can be used for simultaneous spatial and temporal 
beam shaping. This is because, the frequency harmonics 
in our scheme are, at the same time, the angular harmon- 
ics. To start analyzing problem of spatio-temporal beam 
shaping it is instructive to get some analytical informa- 
tion about the phase evolution under the condition of two 
field excitation. First, we assume that E n = f n (r,z)e Ln9 
and that the fields apart from n = 0, 1 are initially zero. 
Then by neglecting the dispersion (J3 n = 0) and the r- 
derivatives of f n (i.e., the radial diffraction), we can find 
an explicit analytical solution to Eqs. (pQ), ([2]). This solu- 
tion gives us z-dependence for all / n 's in terms of Bessel 
functions [10[ . Using known identities for the Bessel func- 
tions one can find the expression for the total field 

x h{z = 0) + f (z = 0)e- i <t- Kz+e A , K = m/c. 

where % = q(z = 0). In our case, the known result for 
E to t [10] is altered by the presence of inside both the 
phase chirp and amplitude modulation. From the am- 
plitude modulation term one can see that the emerging 
intensity pattern of \E to t\ corresponds to a helix. The 
helical shape of the isointensity surfaces can be seen in 
both the (x, y, t) and (x, y, z) subspaces. The time pe- 
riod of the first is 2tt and the spatial period of the second 
is 2tt/K. Changing li to —1 changes the sign in front 
of 0. One can also note, that the directions of rotation 
(chirality) of the helices in t and z are opposite. 

The phase chirp in (|3]) is induced by the Raman co- 
herence and plays a crucial role in the evolution of the 
field, when dispersion and diffraction are taken into ac- 
count. In typical cases the background dispersion of the 
gas is normal, meaning that the compensation of the Ra- 
man chirp occurs for negative detunings from the Raman 
resonance (j = —1). Since diffraction in general and, in 
particular, angular diffraction (~ 8%) is mathematically 
equivalent to anomalous GVD, it seems at the first glance 
that the simultaneous temporal and azimuthal compres- 
sions for materials with normal GVD are not possible, 
because the coherence induced azimuthal and temporal 
chirps have the same sign j. The reality, confirmed by nu- 
merical modeling, is, however, more intriguing. Indeed, 



the angular and frequency harmonics in our case have the 
same complex amplitudes f n . This implies that strong 
coupling between spatial and temporal degrees of free- 
dom occurs in our system. Therefore we can assert that 
if the characteristic length L comi over which maximal 
pulse compression is achieved, is shorter than the char- 
acteristic diffraction length L^, then the temporal and 
azimuthal compression should develop together. Model- 
ing Eqs. $1]) with a plane wave excitation and qo — 0.5 
we found that L com ~ 20 (corresponding to a physical 
distance of ~ 2 — 4cm), which agrees with Refs. 1(1 Il2j. 
In our dimensionless notations is readily estimated as 

To observe the effect of simultaneous temporal and az- 
imuthal compression with j = — 1 (defocusing nonlin- 
earity) and normal GVD (p2 > 0) we take the initial 
excitation like in the modeling shown in Fig. 1, but use 
wider beams (wq = 10\/20, wi = 10-s/lO), in order to 
minimize the role of diffraction over the propagation dis- 
tances ~ L com . The results of this modeling are shown 
in Figs. 2,3. One can see that for a propagation dis- 
tance of z = 20 there are around 30 harmonics gen- 
erated, see Fig. 3(a), implying the presence of vortex 
beams with \l n \ = 15. Simultaneously, a high degree of 
the temporal and azimuthal compression is achieved, cf. 
Fig. 2 (z=20) and Fig. 3(b). With further propagation, 
see Fig. 2, the pulse starts to spread out again produc- 
ing complicated spatial patterns. However, the compres- 
sion in this scheme is a quasi-periodic process, and for 
z = 97 we observe a clear signature of the second, less 
pronounced, temporal and azimuthal compression. The 
spatio-temporal (x,y,i) plot in Fig. 3(c) shows helical 
structure of the intensity \E tot \ 2 for the propagation dis- 
tances before strong defragmentation of the transverse 
profile begins. If we change the detuning sign to j = 1 
(focusing nonlinearity ) , then the spatiotemporal dynam- 
ics is qualitatively different. The compression effect is 
absent and with the same initial conditions we observed 
small scale spatial modulational instability and subse- 
quent self-focusing of the emerging filaments developing 
at z > 30. 

Self-focusing nonlinearity balanced by diffraction sug- 
gests the possibility of spatial solitons. Indeed, the lat- 
ter have been recently reported for the simplest case 
of a two-component Raman model [Hj], while Raman 
self- focusing has been seen experimentally [lf|. The 
question, which is relevant for our problem, is whether 
multi-frequency solitons carrying OAM can be found. 
To find these structures we assumed that E n (x,y,z) — 
f n ( r ) e ine+i(Ki+nK 2 )z^ substituted it into Eqs. Q and 
solved the nonlinear system of ordinary differential equa- 
tions for f n (r) numerically. are the soliton pa- 
rameters chosen to ensure decay of the soliton tails for 
r — > oo. The boundary conditions used at r = are 
lim r ^o /n^~' n ' = const n . Using this approach we have 
found a two-parameter (^1,^2) family of solitons, one 
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FIG. 4: (a) Transverse profiles of \E n (x,y = 13) | 2 for the 
multi- frequency soliton with /i = 1, j = 1, p2 = 0.005. Soli- 
ton parameters are: n\ = —0.5, ^2 = 0. (b) Corresponding 
evolution of the total intensity \Etot\ 2 along z with 1% noise 
in the initial conditions. Isointensity plot at 80% of maxi- 
mum, (c) The same as (b), but for /x = 5. (d) Isointensity 
plot at 20% of maximum for the case of the initial excitation 
with two beams, j = 1 and the other parameters as in Fig. 
1. In (b),(c) and (d) the total intensity is evaluated for t = 0. 



been included (incoherent interaction). In this case, the 
OAM of individual components can be arbitrary, i.e. it 
is not controlled by any selection rules, and higher order 
vortices can not be generated by the system itself. 

In summary: We have demonstrated that a cascaded 
four- wave mixing process in an off-resonantly excited Ra- 
man medium with one of the two input fields being an 
optical vortex leads to generation of multiply charged op- 
tical vortices. Each newly generated vortex beam has its 
own frequency creating strong dependence between the 
spatial and temporal degrees of freedom, which can be 
used for new forms of optical wave synthesis. In particu- 
lar, we have demonstrated the generation of simultaneous 
azimuthally and temporally compressed pulses on the de- 
focusing side of the resonance and helical optical solitons 
on the focusing side. The suggested technique and ob- 
served effects pave the way for practical implementations 
of spatial wave shaping using methods of spectral control. 
The above concepts are also likely to be applicable for the 
waves of a different nature. 



example of which is shown in Fig. 4(a). The total in- 
tensity corresponding to this solution has a helical pro- 
file in both (x, y, t) (not shown, but it is similar to Fig. 
3(c)) and (x,y,z) subspaces, see Fig. 4(b). Solitons for 
sufficiently small values of \i (implying sufficiently large 
values of the intensity) are robust with respect to small 
perturbations, see Fig. 4(b). Taking larger values of \i 
leads to excitation of internal oscillations and wobbling 
helical trajectories, see Fig. 4(c). Notably, experimen- 
tally realistic excitation, by just two beams with Iq = 0, 
Zi = 1, leads to the emergence of structures very close to 
the found soliton solutions, see Fig. 4(d). This is provid- 
ing that the input beams are sufficiently narrow, so that 
the diffraction length is short in comparison to the dis- 
persion length and the self-focusing kicks in before GVD 
desynchronizes the phases of the harmonics. 

The helical soliton beams reported here are qualita- 
tively different from the so-called spiralling solitons or 
rotating soliton clusters [17[, which sustain their rota- 
tion due to the interaction between the individual beams 
accompanied by the conservation of the angular momen- 
tum. In our case the helical evolution does not require 
the presence of more than one intensity lobe and orig- 
inates from the interaction of multiple frequency har- 
monics carrying progressively growing OAM. Ref. [l8[ 
studied complex spatial patterns emerging from the lin- 
early interacting vortex beams carrying different OAM, 
but having identical frequencies. Multicomponent spa- 
tial solitons carrying OAM have been reported in [19[ 
for a nonlinearity which does not depend on the relative 
phases of the interacting harmonics, so that the four- 
wave mixing mediated interaction of the beams has not 
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